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SUBSONIC FINITE ELEMENTS FOR WING-BODY COMBINATIONS 

James L. Thomas 
NASA Langley Research Center 

SUMMARY 



Capabilities, limitations, and applications of various theories for the 
prediction of wing-body aerodynamics are reviewed. The methods range from 
approximate planar representations applicable in preliminary design to surface 
singularity approaches applicable in the later stages of detail design. The 
available methods for three-dimensional configurations are limited as inviscid 
solutions with viscous effects included on an empirical or strip basis. 


INTRODUCTION 


Current research efforts directed toward the design of fuel-efficient air- 
craft dictate that adequate tools be available for the assessment of aerodynamic 
loads across the expected speed envelope. Ashley and Rodden (ref. 1) have sum- 
marized the available methods for aerodynamic analyses of wings and bodies in 
steady and oscillatory motion at both subsonic and supersonic speeds. The ana- 
lytical methods applicable to generalized configurations vary over a range of 
sophistication, accuracy, and computer times required but are generally limited 
as inviscid solutions. Some inviscid-viscid coupling techniques in two dimen- 
sions have yielded good results (refs. 2 and 3), and their inclusion on a strip 
basis into three-dimensional inviscid solutions may serve as a near-term solu- 
tion. The inclusion of viscous effects for generalized configurations across 
the Mach number range remains a far-term solution requiring extensive computer 
resources and advances in turbulence modeling (ref. 4). Immediate design and 
verification methods are thus a combination of experimental and analytical tech- 
niques. The analytical methods largely remain inviscid solutions guided by the 
inclusion of viscous effects on a semiempirical or scrip basis. 

The purpose of this paper is to summarize the capabilities and limitations 
of the existing methods for the steady subsonic analysis of wing-body combina- 
tions. Solutions to the linearized perturbation potential equation (Laplace’s 
equation), with Mach number effects included by the Prandtl-Glauert transforma- 
tion, are considered. Since the governing partial differential equation is 
linear, the solutions may be approximated by distributing a finite number of 
elemental solutions over the body and solving for their relative strengths by 
imposing proper boundary conditions; for example, the flow field must satisfy 
the tangential requirement on the body surface and the Kutta condition at sub- 
sonic trailing edges. Such finite-element solutions have proven to be most 
useful and versatile at subsonic as well as supersonic speeds. The quality of 
the resulting solution is, however, a function of the type, distribution, and 
number of elemental solutions assumed. They require considerably less computer 
resources than the equivalent three-dimensional finite-difference solutions 
required at transonic speeds where the governing equations are nonlinear 
(ref. 5). 
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SYMBOLS 


A 

b 

C L 

C 

P 

c 

c 


aspect ratio, b^/S 
wing span 
lift coefficient 
pressure coefficient 
chord 

section lift coefficient 


1 

d body diameter 

L body length 

M Mach number 

r body radius 

S wing area 

X,Y,Z axis system 

x,z distances along X- and Z-axes 

a angle of attack 

H distance along semispan 

A sweep angle 

A taper ratio 

Subscripts : 

av average 

max maximum 

00 free stream 

f fuselage 


GENERAL SLENDER BODY AND PLANAR WING SOLUTION 

A large number of methods exist for the analysis of planar lifting surfaces 
which account approximately for the presence of bodies. Generally, the methods 
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treat the body separately in an initial analysis and then modify the analysis 
of the lifting surface such that the normal wash on the wing from the body is 
included and the flow is diverted around the body. 

Slender body theory is used in the initial analysis of the body since its 
accuracy is consistent with the assumptions to be made in the wing-body inter- 
actions. Slender body theory assumes the total potential can be composed of a 
far-field potential dependent only on the area distribution and the Mach number 
and a near-field constant-density cross-flow potential solved subject to the 
three-dimensional boundary conditions of flow tangency at the surface (refs. 6 
and 7). The equivalence rule extends the formulation to bodies of general cross 
section as indicated in figure 1. The flow around the actual body differs from 
that of the equivalent body of revolution by only a two-dimensional constant- 
density cross-flow potential that satisfies the flow tangency condition at the 
surface. 

The constant-density cross-flow potential can be solved by any two- 
dimensional method. Dillenius, Goodwin, and Nielsen (ref. 8) have developed a 
solution applicable to noncircular fuselages composed of polar harmonic and two- 
dimensional source-sink terms. A conformal transformation and a distributed 
singularity approach are shown in figure 2. The conformal transformation is at; 
adaption of the Theodorsen technique for airfoil design and was developed by 
Bonner of Rockwell International (ref. 9). The actual body is mapped into a 
circle and the potential for a source or doublet satisfying the boundary condi- 
tions for the equivalent body is transformed back to the physical plane. The 
method is very fast and simple but is limited to bodies in uniform flow fields 
that can be described in polar coordinates as a single-valued function of radius 
versus subtended angle. The distributed singularities approach was developed by 
J. Werner and A. R. Krenkel of Polytechnic Institute of New York and solves for 
the strengths of constant-strength source segments around the body by satisfying 
the flow tangency requirement. The method is applicable to very arbitrary 
bodies in nonuniform flow fields. Comparison of the conformal transformation 
technique of Bonner with experiment (ref. 10) for a parabolic body of revolution 
of fineness ratio 12 and elliptic cross section is shown in figure 3. The 
agreement at this high subsonic Mach number at angles of attack of 0° and 4° is 
generally very good. 

Giesing, Kalman, and Rodden (ref. 11) and Dillenius, Goodwin, and Nielsen 
( *er. 8) have developed methods based on general slender body theory in comb ina- 
xon with vortex-lattice theory and the method of images. In both methods, the 
influt'^e of the body on the Mfting surface is accounted for by including the 
normal wash exterior to the body and then imaging the external singularities 
inside the body. Since the method of images is based on a two-dimensional ana- 
lysis, it does not entirely negate the normal wash from the wing onto the body. 
Thus, the body loading in the nonuniform flow field of the lifting surface and 
imag* system must be recalculated to solve for this residual potential. The 
complete solution is an iterative process in which the continued interaction 
between the body and the lifting surface needs to be computed. However, refer- 
ence 8 has indicated the method is strongly convergent and most of the effects 
are included after the first iteration. The method of images is very attractive 
in that no new unknowns are introduced into the solutions since the •'mage 
strength and location are directly related to the external singularity strengths 
and the geometry of the body cross section. 


i 


i 


i 


13 



An alternate approach has been used by Spangler, Mendenhall, and Dillenius 
(ref. 12) and Woodward (ref. 13) to approximately account for interference 
effects. In their analysis, interference panels are placed on constant-section 
stream tubes of the body. The normal wash from the body is included on the 
lifting surface exterior to the body and the interference panels exist to cancel 
the normal wash induced on the surface of the body. The net result is exactly 
the same as that using the method of images in that the initial influence of the 
body on the wing is Included and the normal wash onto the body from the wing is 
negated. However, there are more equations to solve when the interference 
panels are used, although the region of influence of the wing on the body can 
generally be assumed to be within a couple of chord lengths of the wing root. 

A schematic of the utilization of general slender body theory with a traditional 
vortex-lattice system is shown in figure 4. 

The methods of images (ref. 14) and interference paneling in combination 
with a vortex lattice are compared with an earlier modified Multhopp lifting- 
line approach (ref. 15) for a high-aspect-ratio wing-body combination in fig- 
ure 5. Both the method of images and the method of interference panels give 
similar results and give lower results for the loadings than the earlier 
Multhopp results. Reference 11 has compared the method of images with the 
interference paneling used by Woodward and the agreement is excellent. 

The assumption with either approach is that the flow field around the body 
in the presence of the wing is the same at. that for the body alone. Thin-wing 
assumptions are used which do not account for the finite regions of intersection 
between a wing and a body or the longitudinal acceleration of flow over the body 
on the wing. Because of the singularities trailing downstream with either 
images or interference paneling in accounting for interference effects, the body 
representation is restricted to constant-section cylinders. The methods thus 
give identical results for equivalent positions of the wing above or below the 
midwing position as indicated in figure 6. The results presented are for a 
high-aspect-ratio wing-body combination using a vortex lattice with interference 
paneling. 

A comparison of the theoretical and experimental (ref. 16) span loads for 
a wing-body combination is given in figure 7. All the theoretical methods over- 
estimate the span loading because of the low Reynolds number of the experiment 
(0.3 x 10^). The more approximate theories, however, agree well in the loading 
prediction with the more exact surface singularity representations, such as 
those of Labrujere (ref. 17) or Hess (ref. 18), and, in general, adequate pre- 
dictions of lift and moment are possible with the approximate theories. 

The assumptions of the methods which limit their applicability to general- 
ized configurations also enhance their capability as a preliminary design tool. 
Most of the wing-body interactions are handled and the computer resources 
required are small because of the relatively small number of unknowns. Since 
planar representations are used, the intersection of the wing and body is a line 
and the geometry can be input rapidly. The capability is provided to predict 
quickly and accurately overall lift, moment, and induced drag for complete con- 
figurations at the early design stage, such as in the store separation studies 
of reference 8. The prediction of optimum trimmed loadings subject to lift and 
moment constraints are also possible from a far-field equivalent-horseshoe- 
vortex Trefftz plane analysis such as in references 19 and 20. 
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QUADRILATERAL VORTEX AND SOURCE LATTICE SOLUTION 


A method which computes the interfering flow fields of both wing end body 
simultaneously while still retaining the linearized boundary condition is that 
of Tullnius (ref. 21). The method distributes a series of constant-strength 
quadrilateral vortices over the surface of the body and in the region of the 
wing near the wing-body intersection region as shown in figure 8. Horseshoe 
vortices are used in regions of the wing away from the wing-body Intersection 
region. A source lattice is distributed over the surface of the wing at the 
quarter-chord and three-quarter-chord of each panel, and the source strengths 
are defined as the local slopes of the thickness distribution independent of 
the wing lift. The influence of the quadrilateral vortex dies off rapidly at 
points away from the quadrilateral because of the canceling effects of adjacent 
sides. Hence, the panels can be extended over the fore and aft regions of the 
body. The analysis has been extended to predict thick wing and pylon-fuselage- 
fanpod-nacelle characteristics at subsonic speeds by placing the vortices along 
the mean camber line of the wing (ref. 22). 

Results of the Tulinius wing-body program are compared with experiment in 
figures 9, 10, and 11 for a swept wing-body combination at a Mach number of 0.60 
and an angle of attack of 4°. The unit span load c c and the longitudi- 
nal distribution of fuselage lift c f ^/^niax are P re< ^^ cte ^ very well by the 

theory (fig. 9) ; the fuselage lift increases rapidly in the region of the wing 
root. The pressure coefficients on the wing at two spanwise stations in fig- 
ure 10 and the pressure coefficients on the body at longitudinal stations just 
above and below the wing in figure 11 are also predicted well. The body pres- 
sures are influenced by the wing primarily in the wing root region, and the 
pressures over the aft end of the body are not predicted because of viscous and 
separation effects. The agreement with theory is expected since the wing is 
relatively thin and attached in the midwing position. 

The method cannot account for the longitudinal acceleration of flow over 
the body on the wing (speed bump effect) or equivalent high and low positions 
of the wing because of the linearized planar boundary conditions. The pressure 
coefficients and not just loadings are predicted so that streamlining and con- 
touring of adjacent surfaces at high subsonic Mach numbers can be accomplished. 
Regions of intersecting surfaces are lines so that geometry description is rela- 
tively easy. The number of equations to solve for the simultaneous quadrilateral 
and horseshoe vortex strength increases in comparison with the slender body and 
planar wing analyses but the quality of the aerodynamic solution is higher since 
the body and wing flow fields are solved simultaneously. 

SURFACE SINGULARITY POTENTIAL FLOW 

In order to account for the full potential interactions between the wing 
and body, a surface singularity technique such as that in references 17, 18, 23, 
24, or 25 must be used. In such a method, the singularities are placed on the 
surface of the wing and body such that the tangency and Kutta conditions are 
satisfied. The type of finite-element modeling used for the lifting surfaces 
has been varied, including (1) constant-strength surface source panels with a 
constant-strength vortex sheet on the surface (ref. 18), (2) constant-strength 
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source panels on the surface with interior vortex sheet (ref. 17), or (3) lin- 
early varying source and quadratically varying doublet distributions on curved 
surface panels (ref. 25). Constant-strength source panels have been generally 
used to model the body with the lifting surface carried through the body in order 
to approximately account for the wing carry-throug' ’ift. 

Such a surface singularity approach accounts for the finite Intersection 
region of a wing and body as well as the longitudinal velocity perturbations of 
the body on the wing. However, the method requires a considerable amount of 
geometry specification to panel a complete configuration as shown in figure 12. 
The quality of the resulting aerodynamic solutions are a function of the par- 
ticular finite elements chosen, their placement on the body, and the number 
chosen. Since matrix solutior times are a function of the number of elements 
cubed, the paneling of complete configurations with a minimum of computer time 
while retaining desired accuracy is a difficult task. Recent advances to relieve 
the dependence of the resulting solution on the aerodynamic paneling chosen and 
to reduce the number of unknowns required have been made in references 25 to 27. 

Results for the Hess surface singularity approach (ref. 18) are presented 
in figure 13 for the A * 6 untapered unswept wing attached in intermediate, 
high, and low positions to an infinite circular-cylinder body - the case con- 
sidered earlier with the approximate theory. The local span loading and total 
lift vary with the relative placement of the wing on the body; the body loads 
are shown as average values since the available version of the computer program 
only outputs pressures and integrated loads for the body. The intersection of 
the wing section with the curved body is another curved region that tends to 
accelerate the flow under the wing in a high wing position and above the wing 
in a low wing position. Since the singularities are on the surface, the local 
velocity increase on the lower surface of the high wing decreases the local 
loading and vice versa. Thus, the surface singularity approach yields differ- 
ences in potential theory for high and low wing placement, whereas, the linear- 
ized planar lifting surface theories do not. However, the integrated values of 
lift differ very little with wing placement, indicating again that the approxi- 
mate theories are able to give reasonable estimates of the total forces and 
moments. 

The surface singularity approach is a detail design tool applicable in the 
later stages of design after the initial planform sizes and locations have been 
determined, such as in the design of cruise overwing nacelle configurations in 
reference 28. The inverse design for the surface singularity approach has been 
completed in reference 24, but the procedure for generalized configurations is 
necessarily lengthy and difficult. The surface singularity approach allows the 
calculation of detailed pressure distributions in regions of adjacent surfaces 
(wing fillets, nacelle-strut intersections, etc.) so that contouring and stream- 
lining for minimum adverse pressure and viscous drag can be accomplished. 

CONCLUDING REMARKS 


Various approximate methods utilizing some variation of general slender 
body theory in combination with a planar lifting-surface representation, such 
as the vortex-lattice method or the constant-pressure panel of Woodward, are 
adequate to estimate the loads, moments, and pressures in preliminary design 
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applications. Such methods require limited computer resources and simple 
geometry input specifications and are well suited to inverse design procedures 
since the number of unknowns are small and the planar boundary conditions are 
retained. The methods are most applicable to midwing cases with constant-section 
cylindrical bodies. 

An extension of the vorcex-lattice method to include a quadrilateral /ortex 
representation of the body solves for the wing and body loads simultaneously. 

No restrictions on body shape or wing shape in the intersection regions are made 
although the thin-wing representation is retained. Regions of intersecting sur- 
faces are curved lines and the geometry input remains relatively simple. With 
the method, pressures in regions of adjacent surfaces are predicted to allow 
contouring and streamlining. The method is also well suited to inverse design 
procedures for the wing in the presence of the body since the camber and thick- 
ness solutions are separate. 

In order to accurately predict the correct potential flow pressures in 
areas of intersecting wings and bodies, a surface singularity approach is needed. 
The surface singularity approach removes all thin-wing and linearized-boundary- 
condition assumptions but more than doubles the number of unknowns to be solved 
and the geometry definition required. The detail pressure distributions in 
regions of intersecting surfaces are available so that adverse viscous effects 
can be minimized. 

Viscous effects are not predicted in any of the methods. For the present, 
empirical or strip analyses must be used, such as in the prediction of viscous 
effects using an infinite yawed-wing analogy in two-dimensional strips along a 
swept wing. The usefulness of all the wing-body theories depend on how well the 
theoretical loadings or pressures can be related to the actual physical situa- 
tion. The nonlinear and viscous effects, such as vortex formation near the 
wing-body juncture or separated flow at higher angles of attack, remains untract- 
able computationally. The viscous calculation for generalized configurations 
across the Mach number range remains a far-term solution. 
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Figure 1.- General slender body theory. 
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Figure 2.- Methods for solving two-dimensional cross-flow potential for 

arbitrary cross sections. 





Figure 3.- Comparison of slender body theory and experiment. 
L/d - 12; = 0.9; a/b - 3. 
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Figure 4.- Schematic of utilization of general slender body in 

vortex-lattice theory. 
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Figure 5.- Spanwise load calculation for wing-body combination. 
M ro *0; A - 5; d/c * 0.72; A = 0; A = 1. 
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Figure 6.- Variation of wing placement for A = 6 wing and infinite 
circular-cylinder combination using vortex lattice with interfer- 
ence paneling. ■ 0; d/c - 1; A * 0; A * 1. 




I 


O EXPERIMENTS. 16) 

— SURFACE SINGULARITY^. 18) 
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Figure 7.- Comparison of theory and experiment for A = 6 wing and infinite 
circular-cylinder combination. RAE 101 section; M = 0; d/c * 1; 
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Figure 8.- Quadrilateral vortex and source lattice analysis and design. 







Figure 10.- Quadrilateral vortex and source lattice wing-surface 
pressure coefficients. * 0.6; a ■ 4°. 















